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Abstract 

Following Marleau jjj , we study an extended version of the Skyrme model to which a sixth 
order term has been added to the Lagrangian and we analyse some of its classical properties. 
We compute the multi-Skyrmion solutions numerically for up to B = 5 and show that they 
have the same symmetries as the usual Skyrmion solutions. We use the rational map ansatz 
introduced by Houghton et al. || to evaluate the energy and the radius for multi-skyrmion 
solutions of up to B = 6 for both the SU(2) and SU(3) models and compare these results 
to the ones obtained numerically. We show that the rational map ansatz works as well for 
the generalised model as for the pure Skyrme model. 



1 Introduction 

Recent mathematical developments within the area of non-perturbative methods have 
established the Skyrme model as the strongest candidate for an effective low energy the- 
ory of quantum chromodynamics (QCD). The model was originally proposed by T.H.R. 
Skyrme || to describe hadron interactions. However, it was mainly ignored, until it was 
shown j|, |5], |H that in the large N c limit, where N c is the number of colours, this non-linear 
theory can describe the low energy limit of QCD. This revived the Skyrme model and 
since then significant progress has been made towards the understanding of its properties 
resulting to a relatively successful description of nuclear interactions. 

The Skyrme model is described by an SU(N) valued field U(x, t) which must satisfy 
the boundary condition U — > / as \x\ — > oo, where / is the unit matrix. This condition 
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ensures finiteness of the energy for any field configuration and it also implies that the 
three dimensional Euclidean space on which the model is defined can be compactified 
into S* 3 . As a result, the Skyrme field U corresponds to mappings from S 3 into SU(N). 
Skyrme's idea was to interpret the winding number associated with these topologically 
non trivial mappings as the baryon charge. 

The model is described by the Lagrangian 

C-sk = ^TrR^ + ^Tr[R^ R v ] [R v , R"\ (1) 

where R^ = idfJJ)U is the right chiral current, F n = 189Mev is the pion decay constant 
and a is a dimensionless parameter. The first term in (|IJ) is the non-linear cx-model and 
one can easily show using a scaling argument that with this term alone static solutions 
cannot exist. The same argument shows that one must add to the Lagrangian terms 
involving higher derivatives. This argument led Skyrme to add the second term, usually 
referred to as the Skyrme term, in ([[]) which is the simplest one that preserves the SU(N) 
and Lorentz invariances. 

The Skyrme model can be generalised by adding terms involving higher order deriva- 
tives in the Lagrangian ([!]) |7|, [l], ||, ||. Doing this, one introduces extra parameters that 
can be tuned in to increase the quality of the Skyrme model as an effective low energy 
limit of QCD. For example in J7|, [K| the sixth-order term was used to take into account 
the cu-meson interactions when computing the central Nucleon-Nucleon potential. In a 
different context Marleau studied the model where a large number of higher order terms 
were included in the Lagrangian |], H |9j and where, to avoid the introduction of a large 
number of extra parameters, the coefficients of these extra terms were all related to the 
coefficient of the Skyrme model. 

In this paper we will consider the simplest possible extension of the Skyrme model i. e. 
defined by the Lagrangian (p]) to which we add the sixth-order term 

£ 6 = c & Tr[R^R v }[R u ,R x }[R x ,R^}. (2) 

The unknown coefficient c§ denotes the strength of this term and will be left as a free 
parameter of the model. This particular choice of a sixth-order term is not accidental 
as it is the only term that preserves the Lorentz invariance and the SU(N) symmetry of 
the model and leads to an equation of motion that does not involve derivatives of order 
higher than two. This is the term that was used in J7j. 

In this paper we will focus our attention on the static solutions of the extended Skyrme 
model and thus consider fields that do not depend on time. It is also convenient to define 
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the dimensionless parameter k = 192cqF 2 ci and to introduce the dimensionless units 



y = xy/2/(aF n )y 1 + a/1 + k so that the energy of the model can be written as 

E = —A J dx 3 ( i TrR 2 + Tr 2 + ^A Tr[R u R 3 ] [R v R k ] [R k , rA (3) 



where A = F n / (4y/2a) y 1 + y/1 + k and A = /t/(l + \/l + k) 2 . The parameter A is the 
energy scale of the model. In what follows it will be convenient to use the dimensionless 
energy expressed in the so-called topological units i.e. 

We have chosen this parametrisation of the model so that AG [0, 1] describes the mixing 
between the Skyrme term and the sixth-order term (fj). When A = our model reduces 
to the usual pure Skyrme model while when A = 1 the Skyrme term vanishes and the 
model reduces to what we refer to in what follows as the pure Sk6 model. 

The Euler-Lagrange equations derived from ([3D for the static solutions are given by 
di (lli — -(1 — A) Rj, [Rj, Ri] — —A Rj, [Rj, R k ][Rk, Ri] ^ = 0. (5) 

As mentioned above, an important property of the Skyrme model is that its field cor- 
responds to a mapping from S 3 into SU(N) and as n 3 (SU(N)) = Z each configuration is 
characterised by a an integer which can be obtained explicitly by evaluating the expression 

B = ^l\ dS 3 ^ Tr R i R k) > ( 6 ) 

which following Skyrme's idea is interpreted as the baryon number. Moreover the following 
inequality holds for every configuration 



E > Vl - XB. (7) 

Our extended Skyrme model depends on three parameters: F n , a and Cq or using the 
dimensionless units, A, k and A. To determine the physical values for these parameters, 
we can evaluate different quantities. As our analysis will be purely classical, we will use 
for this purpose the total energy (^) and the isoscalar mean square matter radius given 
by HI 



p2 _ . 2 . f£°drr 2 p B (r) 
n = < r >i=o- 



J™dr p B (r 
where 



p B (r) = Anr 2 B°(r). (9) 
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Notice that after performing the scaling x — > xy/2/ (aF n )\J 1 + \/l + k we can define the 
matter radius evaluated in dimensionless units as 

R = 1 aF n R. (10) 

One can see from the definition of E and R that the ratio of energy or matter radius 
for different solutions only depends on A. In the following section we will evaluate the 
energy and radius of multi-Skyrmion solutions for the general model and evaluate these 
two quantities with the corresponding value for the single Skyrmion and compare them 
directly to the experimental ratio: 

E B = E B {\) Rb = Rb(X) (n) 

E B =i E B=1 (X) R B =i R B=1 (\) ' 

So far all the studies of the classical properties of generalised Skyrme models have 
been focusing on the properties of the single skyrmion (B = 1) |TI], [7|, [II, M. In section 
2 we compute numerically multi-skyrmion configurations for B = 2 to 5 and compare the 
energy and the radius of these solutions with the experimental values. 

It was shown recently^], |l2|] that multi-skyrmion configurations, ie B > 2, can be 
studied systematically using as an approximation the so-called harmonic map ansatz. In 
section 3 we approximate the multi-Skyrmion solution both for the 577(2) and 577(3) 
model using this ansatz. We compare the results obtained with the numerical solutions 
and we show that the harmonic map ansatz provides a good approximation for the multi- 
Skyrmion solutions of the extended model as well. 



2 Numerical Solutions 



In this section we investigate the multi-Skyrmion solutions of the extended SU(2) Skyrme 
model by solving the static Euler-Lagrange equation (||) of the model numerically. Com- 
puting the static solutions of such a three-dimensional model is rather difficult and requires 
a large amount of computing power. As one has to be very careful when assessing the 
accuracy of such numerical results, we are giving in the Appendix a discussion of the 
numerical methods that we have used. 

To compute the solution numerically, it is more convenient to describe the 577(2) fields 
using a four-component vector of unit length, |0| 2 = 1, which is related to the unitary 
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field by U = <fto I + ir ■ <p where / is the unit matrix and f are the Pauli matrices. The 
expression for the energy (Hj) then becomes 



E 



12tt 2 Jr iYM 
A r 
+ 6 



+ 



1 - A 



- 3 |0 M | 2 (0, ■ K ) 2 + 2 (0 M • <f, v ){<f> K ■ 0^(0, • K )j , (12) 

and the Euler-Lagrange equations derived from (|T3|), after adding a Lagrange multiplier 
to impose the constraint |0| 2 = 1, are given by 



1 + (1 - A) |0,| 2 + -A 



121 - 12 -\(<t>„-<t> K ) 2 



+(1 - A) ((0, • </> V(t ) M 



+ 



1 4 



+ A (^V (0„ • 4>vn) I0k| 2 - 0^ (01/ • 0k) (01/ • 0« M ) - 0z/ • 0^k) (4>u ■ 4>k) 

X I X i \ \ A 1 2 i / i i \ | i 1 2 



y i/K yru fn) 



+ 



+ - 



+ (f) u (0 K • W ) (0„ • K ) + y (0 K • M ) (0,, • KAt ) 
1 



vr-3i0 M r(0, 



K ) 2 + 2 (0 M • <p u ) (0 K • M ) (0, ■ K ) =0. (13) 



To compute the B = 1 solution, we use the so-called hedgehog ansatz 

/ sin/(r) sin8 sin(<£>) \ 

, _ sin/(r) sin^ cos(y?) 

^ ~~ sin f(r) cos 6 

V cos/(r) / 



(14) 



where r, # and tp are the usual spherical coordinates. Plugging ( |T4"|) into ( |I~2"D one minimises 
the energy for the profile function f(r) which then has to satisfy an ordinary differential 
equation. This is a very special case of the harmonic map ansatz discussed in the next 
section, so we will just say at this stage that the solutions are radially symmetric and that 
the A dependence of the energy and the radius of the solutions are given on Figure 1. The 
fact that the energy decreases with A is entirely due to our choice of parametrisation; the 
real quantities one has to look at are the energy and radius ratio (|TT|). 



As described in the Appendix, solving (|T3| ) accurately is rather difficult. For this 
reason the case B = 2 was solved differently It is indeed well known that the usual B = 2 
static solution is axially symmetric [p~3] , [16], [14], |15[ and we found that this is also true for 



the extended Skyrme model. Knowing this, we can reduce the system of equations for 

these solutions to a two-dimensional system by using the ansatz 

/ sin / sin g sin(2<£>) \ 
sin / sin g cos(2<^) 
^ sin / cos g 

V cos/ / 



(15) 
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Figure 1: E(X) and R(X) for the B = 1 solutions. 

where = atan(y/x). The two profile functions, f(p,z) and g{p,z) are functions of 
the usual axial coordinates p = \/x 2 + y 2 and z and they satisfy the following boundary 
conditions: 



/(0,0) = 7T f(p -»• oo,^ -> oo) = / P (0,z) =0 

(7(0, z<0) =0 g{0,z>0) = 7T ^U^oc =0 



where i? 2 = r 2 + z 2 . 



(16) 



Substituting (fjJD into (pD we get 



E = — 

6n 



+(1-A) 



(/* + f z ) + sin 2 / (^ 2 + gl) + - sin 2 / sin 2 g 

1 sin 2 / S m 2 s[/ 2 + / 2 + sin 2 f(g 2 p + g 2 z )} + sin 2 / (f p9z - f z9p f 
4 



+A 



P 



2 sin 4 / sin 2 g(f p g z - / 2 # p ) 2 



(17) 



and the corresponding Euler-Lagrange equations are given by 



(f PP + fzz + ~/J - ^ sin 2f sin 2 g - ~ sin 2/ (# 2 + # 2 ) 

+ ( X _ A )| 4 sin2 / sin2 # (fp9p + fz9z) + - sin 2 / [f p g 2 z - f z g P g z ) + \ sin 2/ (f p g z - f z g p ) z 
[P P ' * 



4 
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+ ~2 sin2 f sin2 9 ( f P p + - ^/p ) + sin 2 / sin2 £ 



4 



^(/p 2 + / 2 2 )-sin 2 /(^ 2 + ^ 



+ sin 2 / \jp P g\ + f zz g 2 p - 2f zp g z g p + f p g z g zp - / z #p P # 2 - + fzg P g P z} 

A |^ sin2 / sin 2 f sin2 51 (/ps 1 * ~ f z9 p) 2 + ~t sin4 / sin2 9(fp P g 2 z + f P g z gzp - 
-2fz P g P gz - fzg PP g z + f zz g 2 p + fzg P g P z - f P gz Z g P - - p {f P g 2 z - f z g P g z ) 
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sin 2/ sin 2 / sin 2 g (f p g z - f z g p f 



and 



f 1 \ sin2/ 

g PP + 9zz + -9 P + -r-5-7 
v P J sin7 



2 f 4 

+ + ^2" sin 2 ^ + ( X - A ) I ^2" sin2 / sin 2 9 (dp + 9z) + 

sin 2 # (2 sin 2/ (f p g p + f z g z ) + sin 2 / {g pp + g zz - ) + ^ (/ 2 # p - + / 2 # pp + / 2 ^ 

--^ sin 2#{/ 2 + fl + sin 2 / (g 2 p + g 2 z )} + f z fz P 9p ~ f PP 9zfz + f P f P z9z ~ 2f p g zp f z - j 
+A I sin 2 / sin 2g (f z g p - / p ^) 2 + sin 2 / sin 2 g [f z f zp g p + f z g pp - fp P f z g z - 2f p f z g zp + 



fpfpzgz + f p 9zz ~ fzzf P 9p - ~ (f z 9 P ~ f P fz9z) ~ sin 2 / sin2# (f p g z - f z g p ) z \ = 0. (19) 



P 



2 



2 f „:„ n„ ( f „ t „ ^ 



The advantage of having a two-dimensional system is that we can use much larger grids 
and obtain much more accurate results. As discussed in the Appendix, we have also 
compared the B = 2 solutions obtained by solving (|T3|) and (|18D , (|T9"|) in order to evaluate 
the accuracy of the method we used to solve (p~3|) numerically. 

In figures 2 to 5, we present the A dependence of the energy and radius ratio for the 
B = 2 to B = 5 multi-Skyrmion solutions. We see that in each case the energy ratio 
decreases when the coefficient of the sixth-order term increases while on the other hand, 
the radius ratio increases thus making the multi-Skyrmion solution broader in all cases 
except for B = 2. Tables 1 and 2 compare the energy and radius ratio of the pure Skyrme 
and the pure Sk6 models with the experimental values. We notice that the predicted 
values for the energy are smaller than the experimental values and that the addition of 
the sixth-order term makes the energy ratio even smaller. On the other hand, the addition 
of the sixth-order term makes the multi-Skyrmion solution broader, except when B = 2, 
but the actual values are still much smaller than the experimental ones. 

Another observation we made is that the symmetries of the multi-Skyrmion solutions 
for the general model are the same as for the pure Skyrme model. The solutions for 
5 = 2,3 and 4 Skyrmion have respectively the shape of a torus, a tetrahedron and a cube 
while the B = 5 Skyrmion solution has the same D2d symmetry. 

It is a well know problem that the binding energies predicted by the Skyrme model are 
too large and that the radius of the classical solutions are too small. One usually argues 
that quantising the model will somewhat solve this problem. Adding the sixth-order 
term does not improve this: the energy binding is even stronger and the multi-Skyrmion 
solutions are narrower except for B=2. 
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3 Harmonic map ansatz 



In this section we will use the rational map ansatz to compute configurations that ap- 
proximate solutions of the extended Skyrme model. We will then use these configurations 
to evaluate the energy and radius of the multi-Skyrmion configurations, check how these 
properties depend on A and compare these results to the ones obtained for the numerical 
solutions. 

The rational map ansatz, introduced by Houghton et al. || is an extension of the 
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Figure 4: E(X) and R(X) ratio of B = A/B = 1 for the numerical solutions. 
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Figure 5: E(X) and R(X) ratio of B = 5/B = 1 for the numerical solutions. 

hedgehog ansatz found by Skyrme, which using the usual polar coordinates is given by 

U(r, 6, if) = exp(ig(r) n(0, tp) > tr). (20) 

In the hedgehog ansatz n is a unit length vector describing the one-to-one mapping of 
the two-sphere into itself and g(r) is a profile function satisfying the boundary conditions 
c/ (0) = 7r and g(oo) = 0. The rational map ansatz consists in using for n(0,ip) harmonic 
maps from S* 2 into S 2 while keeping the same boundary conditions for g. One can easily 
show that the baryon number for such a configuration is given by the degree of the 
harmonic map. To approximate a solution of a given baryon charge, one takes for n the 
most general rational map of the given degree, inserts the ansatz into the expression for the 
energy and tries to minimise this expression with respect to the parameters of the rational 
map and the profile functions g. When doing so, the integration over the radius r and the 
angular variables 9 and (p decouple and the rational map appears only in two expressions 
integrated over the whole sphere. One of them can be evaluated explicitly and is equal to 
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Experiment 


Numerical solutions 


R 

1 J 


Enerev (Me VI 


Ratio 


Skyrme Ratio 


Sk6 Ratio 


1 


939 








2 


1876.1 


1.99798 


1.9009 


1.8395 


3 


2809.374 


2.99188 


2.7650 


2.7103 


4 


3728.35 


3.97055 


3.6090 


3.5045 


5 


4668.795 


4.97209 


4.5000 


4.3780 



Table 1: Experimental energy ratio, Eb/Eb=i, and values obtained for the numerical solutions. 



The experimental values (MeV) correspond to isotopes with minimum mass [17]. 





Experiment 


Numerical solutions 


B 


Radius (fm) 


Ratio 


Skyrme Ratio 


Sk6 Ratio 


1 


0.72 








2 


1.9715 


2.73819 


1.3549 


1.308 


3 


1.59 


2.2083 


1.5080 


1.5570 


4 


1.49 


2.06944 


1.6850 


1.7420 


5 






1.8890 


1.9250 



Table 2: Experimental radius ratio, Rb/Rb=i and values obtained for the numerical solutions. 
The experimental values (fm) correspond to nuclei with minimum mass. [11, 18, [H], 2C] 

the degree of the harmonic map while the other must be minimised with respect to the 
parameters of rational map. Doing so leads to a unique rational map, up to an arbitrary 
rotation, which describes the radial dependence of the Skyrmion configuration. Then 
one minimises the effective energy by solving the Euler-Lagrange equation for the profile 
function g. The configurations obtained by this construction have the same symmetries 
as the exact solutions and their energies are only 1 or 2 percent higher |2|. 



This construction was later generalised by Ioannidou et al. [12 to approximate so- 
lutions of the SU(N) Skyrme model using harmonic maps from S 2 into CP N ~ l . The 
generalised ansatz takes the form 



U(r,9,<p) 



,2ig(r)(P(B, V )-I/N) 



-2ig(r)/N 



(/ + (e 2i9{r) - l)P{ 



(21) 



where P(6, <p) is an N x N projector. As we want to study some solutions of the SU(3) 
model as well, we will use the generalised construction. At this stage it is convenient to 



10 



introduce the complex coordinate £ = tan(9/2)e l ' p which corresponds to the stereographic 
projection of the unit sphere onto the complex plane. 



The procedure to minimise the energy is the same as the one outlined above where 
the projector will be taken as a harmonic map from S 2 into CP N ~ X i.e. a projector of 
the form pTJ 

m 

where / is a N components complex vector whose entries are all rational functions of £. 
The degree of the harmonic map is given by the highest degree of the components of / 
and the baryon number is again given by that degree. 



P(f) = ^ ( 22 ) 



Substituting the ansatz (|2~ID in the general expression of the energy density (f|) we get 



E = — dr (A N g 2 r r 2 + 2M sin 2 g (1 + (1 - \)g 2 ) + (1 - \)X- 
on J 



r z 6 r 4 



where : 



An = ^(N-l), 

M = ^|^^(l + |£| 2 ) 4 Tr([^P,%P] 3 ). (24) 

The integral M is nothing but the energy of the two-dimensional Euclidean CP N ~ l 
ex-model and for the harmonic projector it is equal to the degree of the harmonic map, 

B=^-jdt di Tr (P [%P, 9 e P]) . (25) 

As X and M. are independent of r they can be minimised with respect to the parameters of 
the harmonic maps. In what follows we will prove that M. is identically zero so only X will 
have to be minimised, something which was already done in 0] and ||12|| . The minimisation 



of the energy with respect to the profile function g{r) is then straightforward. 

To prove that M. vanishes, we need to use some properties of the projectors given by 
(f22|) where df/d£ = 0. First of all, it is easy to check that 

PP f = and P^P = P^ (26) 
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where P^ denotes the derivative of P with respect to £ and from this we have 

PP^P = and thus Pf = 0. 



(27) 



Using (|27| ) we notice that 



Tr[Pt,Pii 



Tr(P^ - P^) 



Tr {(P^r + (-l)"(^n 
(l-(-l) n )Tr ((P f J£ 



proving that 



(2f 



(29) 



As a result, the energy density ( ]2"3| ) simplifies further and if we treat M and X as two 
parameters then one can minimise the energy E by solving the following Euler-Lagrange 
equations for g 

(-. . r, K r l - X sin2 9 . T A sin 4 #\ 2 / A sin 4 ^\ 

p rr l + 2jv— ^+J- —)+-g r \l-l- 



Tr[P^,P^\ n = for n odd, 
and thus that M. in ([24] ) is identically zero. 



1 sin 2t7 



A N r 4 



A N r 4 



sm g 



An 



tf((l-\))gi-l)+I—JL(\gi-l + \) =0. (30) 



We see from our analysis that the harmonic maps for the extended Skyrme model are 
the same one as the usual Skyrme model. The harmonic map ansatz predicts thus that 
the solutions of the usual and the extended Skyrme models have the same symmetries. 
This has been confirmed by the numerical solutions. The only difference, for the ansatz, 
between the two models comes from the profile function. This is due to the presence of 
the extra terms appearing in ( j3Cf) . 

In Table 1 we list the minimum energy for the harmonic maps that we will use later 
together with the corresponding value for X. 



3.1 Energy and radius ratios for the SU(2) model 

In this section, we analyse how the properties of the multi-Skyrmions rational map ansatz 
depend on the parameter A. Using the value of I given in Table 1, we compute the profile 
g by solving ( j30|) and evaluate both the total energy and the radius of the configurations. 
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Harmonic Map /(£) 
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ft i \t 
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Cs» !) 
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5.81 


/ ^ , — \ i 
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4 


Q 
O 


\ S \S V ^"ji V ^ L s / 


13.58 


ff 3 , 1.576f, v^ -1 )* 


10.65 


4 


+ 2V3«r + 1,4- 2V3< Z + 1 J 


on /? r 

20.65 


(4 , 2.71914^, 1) 


18.05 


5 


(e(e 4 + ^ 2 + a),< 4 -^ 2 + i) 4 , 


37.75 


(e 5 -2.7e,2e 4 + i,9/2a* 


27.26 




a = 3.07, 6 = 3.94 








5* 


(£(£ 4 -5),-5£ 4 + l)* 


52.05 






5** 


(4 5 , 1)* 


84.425 







Table 3: Harmonic maps /(£) minimising the angular integral X for SU(2) and SU(3) 
The 5* and 5** configurations denote saddle points that we also consider. 

In Figures 6 to 9, we show the energy ratio and the radius ratio defined in ( pT|) for different 
values of the baryon number. At this stage we would like to remind the reader that A = 
corresponds to the pure Skyrme model while A = 1 is equivalent to the pure Sk6 model. 
Moreover, the ratio presented on the figures only depends on A i.e. the mixing between 
the two Skyrme terms. When B > 6, the graphs we obtained were all similar to Figure 5. 

When comparing these results with the numerical solutions, we notice first of all 
that the energy ratio predicted by the ansatz is always too large. Apart from this, the 
prediction for the energy is rather good except for the case B = 2 where the energy 
difference between the numerical solution and the rational map ansatz is 7 times as large 
for the Sk6 model than for the pure Skyrme model. 

We also notice that the graphics obtained for the numerical solutions do not exhibit 
local minima as observed on the graphs obtained for the harmonic map ansatz. The only 
exception is the radius ratio obtained for the B = 5 exact solution but the effect is so 
small that it could be a numerical artefact. 

The radius ratio obtained with the harmonic map ansatz is always too large when 
compared with the radius ratio of the exact solution. For B = 2, the radius ratio increases 
with A and the error only gets worse as A increases. The case B = 3 is rather surprising 
as the radius ratio has a deep local minimum around the value A = 0.3; this is where 
the relative error is the smallest, otherwise the relative error is smaller for the pure Sk6 
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model than for the pure Skyrme model. The cases B = 4 and B = 5 are very similar: 
the radius ratios decrease when A increases and the error for the pure Sk6 model is very 
small especially when B = 4. 

We can thus conclude that the harmonic map ansatz produces good approximations 
to the solutions of the generalised Skyrme model and the error is in most cases smaller 
for the pure Sk6 model than for the pure Skyrme model, the only exception being the 
case B = 2. 




Figure 6: E and R ratio of B = 2/B = 1 as a function of A. 
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Figure 7: E and R ratio of B = 3/B = 1 as a function of A. 



So far we have examined the behaviour of the model for harmonic maps that minimise 
the angular integral X and correspond to minimum energy configurations. We have next 
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Figure 9: E and R ratio of B = 5/B = 1 as a function of A. 

considered harmonic maps that correspond to saddle points of the energy for B = 5. 
The reason behind this selection lies on the fact that the binding energies of the multi- 
Skyrmion solutions are much larger than the experimental values. For this reason, we 
have considered the two harmonic maps, B = 5* and B = 5**, given in Table 3. The first 
one has octahedral symmetry whereas the second gives a toroidal Skyrme field || . 

For the case of B = 5*, shown on Figure 10, we see that the binding energy is slightly 
larger than 5 for the pure Skyrme model and that it decreases when the strength of the 
sixth-order term increases, going through the experimental value 4.97 when A ~ 0.1. 

The second case, shown on Figure 11, is the only example where we have seen a local 
maxima for the energy ratio that is larger than the energy of both the pure Skyrme and 
Sk6 model. 
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3.2 Energy and radius ratios for the SU(3) model. 



In this section we look at the harmonic maps configurations for the 577(3) models |T2 . 
The harmonic maps that we will use and the corresponding values of X are all given in 
Table 3. The single ££7(3) skyrmion is the well-known hedgehog ansatz and it is just an 
embedding of the £77(2) solution. 



Notice that the numerical constant appearing in (|30|) and (|23|) is now equal to 4/3. 

We should stress here that these configurations approximate solutions that are believed 
to be saddle points of the energy. Their energy is larger than the corresponding SU(2) 
embeddings and they have a different symmetry as well. 




Figure 12: E and R ratio of B = 2/B = 1 as a function of A for the SU(3) harmonic map 
ansatz. 



It is interesting to notice that unlike the SU(2) model, the energy of the B = 2 
solutions increases with A. For a given B and a fixed value of A, the energy ratio of these 
configurations is always larger than the energy ratio of the corresponding ££7(2) solutions, 
while on the other hand, the radius ratios is always smaller. 

It is also interesting to notice that the A dependence of the energy and radius ratios 
obtained for a given B looks very much like the curve obtained for the SU (2) model for 
B — 1 Skyrmions. This can be explained by performing the change of variable r — > rk, 
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2.96 




E and R ratio of B = 3/B = 1 as a function of A for the SU(3) harmonic map 




Figure 14: 
ansatz. 



E and R ratio of B = 4/B = 1 as a function of A for the SU(3) harmonic map 




Figure 15: E and R ratio of B = 5/B = 1 as a function of A for the SU(3) harmonic map 
ansatz. 
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Figure 16: E and R ratio of B = 6/B = 1 as a function of A for the SU(3) harmonic map 
ansatz. 



where A; 2 = ^(1 + Jl + 4X/(A N (1 - A) 2 )), and rewrite (|| 



as 



/ 4^ /rfr (^^ 2 r 2 + 2Ar'sin 2 ^(l + (l-A')^ 2 ) + (l-A / )X'^# + A'J' SUI " - 2 



3tt 



k 2 r 2 



r 



l 

'/ 

2 SV/ 



(31) 



where M' = N/A^, I 1 = I/An, and A' = A/fc 4 . The function k(X) monotonically 
decreases from k(0) = 1 to k(l) = A^^ A and so it is relatively close to 1 for all values of 
A. In Table 4, we give the values of M' and I' for the SU(3) ansatz and we notice that 
the SU(3) solutions for B = 4 and B = 5 are closely related to the SU(2) solutions for 
respectively B = 3 and B = 4. 



SU{2) 


SU{3) 




N 


1 


N 


1 


Af' 


T 


2 


5.81 


3 


10.65 


2.25 


7.98 


3 


13.58 


4 


18.05 


3 


13.54 


4 


20.65 


5 


27.26 


3.75 


20.44 


5 


37.75 


6 


37.33 


4.5 


28 



Table 4: N 1 and X' for the SU(3) anstaz. 

In tables 5 and 6 we compare the energies and the radius ratios, obtained for the 
SU{2) and SU(3) model using the rational map ansatz. We also compare these values 
with the SU(2) numerical solutions. 
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SU{2) Numerical Solutions 


SU{2) 


517(3) 


D 
D 


SIka/tttip R^tin 


Sk6 Ratio 


kJ -Lv y L lllC -L \iCX u Iv 


Sk6 Ratio 


SikA/TTYIP n^flO 
kj -Lv y x me j. KjCXi u iv 


Sk6 Ratio 


2 


1.9009 


1.8395 


1.96223 


1.95407 


1.98468 


1.99585 


3 


2.7650 


2.7103 


2.88541 


2.82888 


2.95974 


2.94756 


4 


3.6090 


3.5045 


3.69164 


3.52850 


3.84491 


3.76064 


5 


4.5000 


4.3780 


4.65685 


4.45345 


4.72485 


4.56766 


6 






5.54105 


5.26743 


5.57660 


5.33739 



Table 5: Energy ratio, Eb/Eb=i, for the SU(2) numerical solutions, and the SU{2) and 517(3) 
rational map anstaz configuration. 





SU(2) Numerical Solutions 


SU{2) 


SU{3) 


B 


Skyrme Ratio 


Sk6 Ratio 


Skyrme Ratio 


Sk6 Ratio 


Skyrme Ratio 


Sk6 Ratio 


2 


1.3549 


1.308 


1.37023 


1.39403 


1.20691 


1.234384 


3 


1.5080 


1.5570 


1.63107 


1.62894 


1.45842 


1.483996 


4 


1.6850 


1.7420 


1.78911 


1.746286 


1.63002 


1.62755 


5 


1.8890 


1.9250 


2.013822 


1.95551 


1.78149 


1.75505 


6 






2.178298 


2.09768 


1.909141 


1.859916 



Table 6: Radius ratio, Eb/Eb=i, for the SU(2) numerical solutions, and the SU(2) and SU(3) 
rational map anstaz configuration. 

3.3 Conclusions 

We have studied an extension of the Skyrme model defined by adding to the Lagrangian 
a sixth-order term. We have computed the multi-Skyrmion solutions of the extended 
model for up to B = 5 and we have shown that they have the same symmetry as the 
pure Skyrme model. We have analysed the dependence of the energy and radius of the 
classical solution with respect to the coupling constant A. We found that the addition 
of the sixth-order term makes the multi-Skyrmion solution more bound than in the pure 
Skyrme model and that it also reduces the solution radius. 

We have also used the harmonic map ansatz to approximate the numerical solutions 
and we found that the ansatz works as well, and in many cases even better, for the 
extended model than for the pure Skyrme model. 
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A APPENDIX 



In this appendix we describe the numerical methods that we have used to compute our 
numerical solutions. 

A.l 3 Dimensional Solutions 

To compute the three dimensional solutions described in section 2, we discretised the 
static equations using finite differences and we solved them using the relaxation method. 
We used the fixed boundary condition, taking the vacuum value for the field on the edge 
of the grid. 

The values obtained for the energy and the radius with our methods are affected by 
two sources of inaccuracy. The first one is the finiteness of the grid which, by distorting 
the field slightly, increases the value of the energy. The second one is the fact that finite 
differences systematically underestimate the value of the energy. One could of course hope 
that the two effects cancel out, but as it is difficult to evaluate their order of magnitude 
one has to experiment and reduce them both as much as possible. 

To reduce the edge effects, we computed the same solutions on grids of different sizes 
L but keeping the lattice spacing dx = L/N constant, where N is the number of lattice 
points in each direction. We then looked at how the energy changed as a function of the 
size and chose a value for L for which the energy is only slightly affected by the edge 
effects. 

The finite difference scheme we used is of order two, so when we evaluate the total 
energy we can write 

E = E + E 1 dx + E 2 dx 2 + O{dx 3 ). (32) 

In theory, E\ = but in practice it is a small but non-zero coefficient induced by the edge 
effects. To improve the evaluation of the energy for a given solution, we computed the 
solution for at least three different values of dx using a grid of size L for which the edge 
effects are sufficiently small. We then fitted these values to the coefficient Eq, E\ and E 2 in 
(j32|) getting E as a better estimation for the energy. Notice also that E 2 is always negative 
and that |ET|*d:c <C \E2\. When this last condition is not satisfied one must conclude that 
the edge effects are large and one must increase L. To check our evaluation we performed 
the same interpolation for the topological charge Q = Qo + Q\dx + Q 2 dx 2 + 0(dx 3 ). As 
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we know that it must be an integer B, the quantity (B — Qq)/B is a good estimation of 
the relative error on the topological charge but also on the energy Eq. 



For the solution B = 2 ... 5 we used a box ranging from —8 to 8 in all directions and 
we used grids of 100 and 120 and 140 points. We also found that for a given value of B, 
E 2 did not change much with A. We were thus able to evaluate it for a few values of A 
and used an extrapolation for the other values. We also found that the relative error on 
the energy was smaller than 0.5%. 

As an alternative method to evaluate the energy we have considered computing the 
quantity E/Q, as if E 2 and Q 2 were comparable, we would not have to compute the 
solutions for different values of dx. Unfortunately we found that with our discretisation, 
Q 2 is about 50% larger than E 2 and as a result the value we get for E/Q increases when 
dx decreases and thus underestimates the energy value. 

To evaluate the radius, we have used the same method, but for this quantity the 
integrand decreases more slowly towards infinity and as a result the value is more affected 
by the finiteness of the grid. We believe that the overall behaviour of the radius ratio 
graph can be trusted but some of the fine details might be numerical artefacts. 

To double-check our results, we have computed the B = 2 axially symmetric solutions 
by solving (|i~8|) and (|l^) on a two-dimensional grid. This made it possible to use many 



more points and much larger grids. When using the grid defined by z G [—20, 20] and 
r G [0,20] taking dx = 0.05 the error was smaller than 0.1% and we found for example 
E = 2.378 for the pure Skyrme model. When computing this solution by solving the 
three-dimensional equation using the method described above, the difference between the 
two energies was less than 0.1%, thus validating the methods that we used. 

The energy values that we obtained for the pure Skyrme model all fit within 1% the 



value given in [13] except for B = 2 where the error is 1.5%. As the numerical methods 



and the type of grid used are not described in |L3[ it is difficult to make any further 
comparison between the numerical results. 



A. 2 profiles 

To compute the profile functions for the hedgehog ansatz or the harmonic map ansatz we 
have used both the shooting and the relaxation methods. We have in every case compared 
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the solutions obtained with grids of different sizes and different number of points to ensure 
that our results were accurate and that they were not affected by edge effects. We were 
led to use very large grids, R ma x = 80, to get an accurate value for the radius as well as 
up to f60000 lattice points. 
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